Baryonic Strangeness and Related Susceptibilities in QCD 
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The ratios of off-diagonal to diagonal conserved charge susceptibilities e.g., Xbs /xs ,XQS /xs , 
related to the quark flavor susceptibilities, have proven to be discerning probes of the flavor carrying 
degrees of freedom in hot strongly interacting matter. Various constraining relations between the 
different susceptibilities are derived based on the Gell-Mann-Nishijima formula and the assumption 
of isospin symmetry. Using generic models of deconfined matter and results form lattice QCD, it is 
demonstrated that the flavor carrying degrees of freedom at a temperature above 1.5T C are quark- 
like quasiparticles. A new observable related by isospin symmetry to Cbs — —3xbs/xs and equal 
to it in the baryon free regime is identified. This new observable, which is blind to neutral and 
non-strange particles, carries the potential of being measured in relativistic heavy-ion collisions. 

PACS numbers: 12.38.Mh, ll.10.Wx, 25.75.Dw 
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I. INTRODUCTION 

The goal of of the heavy-ion program at the Relativis- 
tic Heavy Ion Collider (RHIC) is the creation and study 
of heated strongly interacting matter at nearly vanishing 
baryon density 0, @ ■ Detailed models of nuclear reac- 
tions predicted that the energy deposition in the center- 
of-mass frame should be sufficient to cause temperatures 
at mid-rapidity in central collisions of gold nuclei to reach 
upwards of 300 McV 3|. These predictions have been 
confirmed by the experimental results of the four RHIC 
detector collaborations, which set a lower bound of about 
5 GeV/fm 3 on the energy density at a time r = 1 fm/c 
in central Au+ Au collisions Q . This estimate for the at- 
tained energy density should place the produced matter 
well into region of the QCD phase diagram that cannot 
be described as a dilute hadronic resonance gas. Until 
recent results from the RHIC experiments have raised 
doubt about this interpretation 0], the matter in this 
domain of the QCD phase diagram was expected to be 
a colored plasma composed of quasiparticle excitations 
with the quantum numbers of quarks and gluons 0, Q . 

Signals of the excited matter produced in the early 
stages of such a collision, buried in the pattern of de- 
tected particles, may be divided into two categories: hard 
probes, such as the modification of partonicjets by the 
medium [(|, bound states of heavy quarks Q, and elec- 
tromagnctically produced particles y|; and bulk obscrv- 
ables, dealing with low momentum particles which make 
up a large fraction of the produced matter. Bulk ob- 
servables include the single inclusive spectra of identified 
hadrons [t| and the event-by-event fluctuations of con- 
served charges 0. The latter are the subject of our 
present study. 

The theoretical analysis of the observed suppression of 
energetic hadron emission in Au+Au collisions at RHIC 
("jet quenching" 0) confirms that matter with a very 
high energy density is produced. This matter clearly ex- 
hibits collective behavior as evidenced by its radial and 
elliptic flow 0]. As a large elliptic flow requires large 
pressure gradients which can only be present during the 



earliest stage of the collision, the matter produced at such 
times must exhibit the properties of a fluid. A quantita- 
tive analysis of the observed magnitude of the flow indi- 
cates that this fluid must be nearly ideal, i. e. endowed 
with a very low viscosity 0|. This result suggests that 
there must exist a strong interaction between the con- 
stituents of the medium. 

At vanishing baryon density, the entire range of ther- 
mal conditions expected to be attained at RHIC may 
be simulated by the numerical methods of lattice QCD 
(LQCD) at finite temperature [I^. Here one com- 
putes the grand canonical partition function of a sys- 
tem whose states are thermally weighted by the QCD 
action at a temperature T, with baryon, electric charge, 
and strangeness chemical potentials set to zero (/is = 
MQ — Ms — 0)- Explorations on the lattice consist of a 
three-fold approach studies of the behavior of the 
components of the stress energy tensor, i. e. the energy 
density e and the pressure P; spatial and temporal cor- 
relation functions, and the recently measured derivatives 
of the free energy such as the various conserved flavor 
susceptibilities. 

Investigations of the first kind, exploring the QCD 
equation of state at vanishing chemical potentials, pro- 
vide the most solid evidence for the expectation that, 
when hadronic matter is heated beyond a critical tem- 
perature T c ~ 170MeV, a transition to a new state of 
matter, the quark-gluon plasma (QGP), will occur. The 
transition is signalled by a steep rise in the energy density 
and the pressure as a function of the temperature. The 
slow rise of both quantities prior to the sudden transition 
has come to be understood in a picture of a hadronic res- 
onance gas 0| . However, attempts at an description of 
the excited phase as a weakly interacting plasma of quasi- 
particles 0, ^| has not met with success in the region 
from T C <T < 3T C . 

This finding might indicate that matter in this region 
may not be a weakly coupled plasma where quarks and 
gluons are deconfined over large distances as it was orig- 
inally proposed ^?J- ^ has been established that such 
a weakly coupled state, indeed, occurs at much higher 
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temperatures (see 01 an d references therein). Assuming 
that temperatures at RHIC do not exceed 3T C , a strongly 
coupled state is not inconsistent with the strong collec- 
tive behavior observed in experiments. It is clear that 
a microscopic understanding of the emergent degrees of 
freedom in this regime is essential for an explanation of 
the "perfect liquid" character of the matter created in 
nuclear collisions at RHIC. 

There already exists a candidate model for such mat- 
ter, proposed by Shuryak and Zahed consisting of 
a tower of colored bound states of heavy quasi-particulatc 
quarks and gluons. In spite of the presence of such heavy 
quasi-particles, this model was to account for the large 
pressure required by the RHIC data because of the prolif- 
eration of excited bound states. The large cross sections 
resulting from resonance scattering were put forward as 
the cause for the short mean free path which is a requi- 
site for the appearance of hydrodynamic behavior. The 
model has, however, fared poorly in a comparison with 
lattice susceptibilities. Koch et al. |2l| recently proposed 
the susceptibilities of conserved charges B (baryon num- 
ber), Q (electric charge) and S (strangeness) and their off 
diagonal analogues as diagnostics of the conserved charge 
or flavor carrying degrees of freedom in a strongly in- 
teracting system. The primary quantity of interest is 
the ratio of the covariance between baryon number and 
strangeness <J 2 BS oc xbs to the variance in strangeness 
a s x XSi renormalized to be unity in a quasi-particle 
gas of quarks and gluons, 
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This quantity may be calculated on the lattice and esti- 
mated in the dynamical model of ref . [l^ ■ The calculated 
values on the lattice were found to be 50% higher than 
those computed in the model |2ll |. Such comparisons 
demonstrate the capacity of ratios such as Cbs to serve 
as tests for models of the QGP. The origin of the dif- 
ference between the results from lattice QCD and those 
from the bound state model will be further discussed in 
the upcoming sections. 

The objective of the present work is to continue the 
study of diagonal and off-diagonal flavor susceptibilities. 
In the next section, we point out that there exists an 
entire gamut of such diagonal and off-diagonal suscep- 
tibilities for conserved quantum numbers, and they are 
related by a simple transformation to the equivalent basis 
of quark flavor susceptibilities. This presentation builds 
on the work of Gavai and Gupta [2^ ■ In the next section, 
a complete set of such susceptibilities is introduced and 
related via the Gell-Mann-Nishijima formula. The alter- 
native set of quark flavor susceptibilities, which appear 
more often in the literature, is also compared. In Section 
III operator relations between the off-diagonal suscep- 
tibilities in different bases are derived and the generic 
behavior of such operators depending on the prevalent 
degrees of freedom is outlined. In Section IV several 
models for the flavor carrying sector of the QGP are ana- 



lyzed, and arguments favoring a picture of quasi-particlc 
quarks are presented. In Section V we formulate observ- 
ables based on the ratios of susceptibilities that may be 
estimated from experimental measurements. Concluding 
discussions are found in Section VI. 



II. DIAGONAL AND OFF-DIAGONAL 
SUSCEPTIBILITIES 



In the lattice formulation of QCD, the fundamental 
degrees of freedom are local quark and gluon fields. Un- 
der conditions where deconfinement has been achieved, 
the elementary set of conserved charges is given by the 
quark flavors: the net "upness" (Am = u — u), "down- 
ness" (Ad = d—d) and "strange-quarkness" (As = s — s). 
An alternate basis is provided by the hadronically defined 
conserved charges of baryon number (B), electric charge 
(Q) and strangeness (S). The two bases are related by 
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In what follows, the "A" will be omitted and the vari- 
ables u, d, s will be understood to denote the net flavor 
contents. 

The mean values of any conserved charge may be mea- 
sured in a thermal ensemble of interacting quarks and 
gluons on the lattice. The grand canonical partition func- 
tion may be defined using either basis i.e., 

Z(T,hb,Hq,Hs) = Z(T, fi u , fi u , n s ), 

V-B - MQ 



Ms = -Ms- (3) 



The mean values and variances of any combination of 
conserved charges may be obtained from appropriate dif- 
ferentiation of either partition function: 



d 

W = T—\ogZ(T, ix x , iiy), 



log Z(T, n x ,Hy). 



(4) 
(5) 



While the mean values in a given basis are related to the 
other via Eq. J2J), the variances exhibit a more complex 
structure. A similar and larger matrix relates the six fluc- 
tuation measures, viz. the variances o b ,Oq,o 2 s and the 
covariances <J 2 BS , ct B q, Cq S , to the six diagonal and off- 
diagonal quantities constructed from the quark flavors. 
The 6x6 matrix relating these two sets of (co-)variances 
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is given by 
( B 2 \ 



BQ 
BS 
QS 

V * 2 J 



1 1 

4 1 

2 -1 





2 2 2 

-4 -4 2 
1 1 -2 
-3 -3 
0-6 3 







ud 
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ds 

\* 2 ) 



(6) 



where the corresponding subscripts of cr 2 y are used to 
indicate the corresponding variance. The above matrix 
immediately demonstrates the utility of using ratios of 
a BS' a QS ano - a s as PP ose d to the other three variances, 
as these form a smaller subgroup with the quark flavor 
covariances a 2 s , o 2 ds and the strangeness variance a 2 . The 
(co-) variances are extensive quantities: 



<*x y = VT Xxy, 



(7) 



where \ xy is the intensive diagonal or off-diagonal sus- 
ceptibility. These susceptibilities can be measured on the 
lattice. In heavy-ion experiments, the variances and co- 
variances are measured by means of an event-by-event 
analysis of the corresponding conserved quantities, i.e., 
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where E represents the set of events, Ne is the number 
of events considered, and Xi , Yi are the net values of the 
conserved charge in a given event i. The volume inde- 
pendent ratios of variances, measured event-by-event in 
heavy- ion collisions, may then be directly compared with 
the lattice estimate for the ratio of susceptibilities. 

Besides relating the susceptibilities or variances from 
one basis to another, simplifying relations may be ob- 
tained between the variances in a given basis using the 
Gell-Mann-Nishijima formula, 



B + S 



(9) 



where, I3 denotes the third component of the isospin. 
Because the mass difference between the u, d quarks is 
small compared with the typical scale of hadron masses, 
the masses of all hadrons in a given isospin multiplet are 
degenerate. One may also assume that any quasi-particle 
excitations in highly excited chromodynamic matter dis- 
play a large degree of isospin symmetry. It should be 
pointed out that in actual experiments, besides a small 
mass split between the u and the d quark, there is an over- 
all isospin asymmetry in the entire system as all heavy 
nuclei have a large isospin due to their proton-neutron 
imbalance. Such an isospin asymmetry is, however, dom- 
inantly correlated with the net baryon number of the nu- 
clei. As the matter formed in the central rapidity region 
has nearly vanishing net baryon density, one may assume 



that the produced matter has almost no net I3 compo- 
nent, on the average. In the following, we will also ignore 
the isospin violating effects of weak decays and electro- 
magnetic interactions. 

To obtain relations between quadratic variables involv- 
ing strangeness, we obtain the variance of the left and 
right hand side of the equation 



QS = hS 



BS + S 2 



(10) 



Taking the ensemble or event average of the above quan- 
tity, we obtain 
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In the above, the total charge and strangeness measured 
in an event i is denoted as Qi, Si. In the case that the 
degrees of freedom or active flavors / are eigenstates of 
charge or strangeness, this expression may be decom- 
posed as Qi — Ylf n iQf wriere n{ are the number of 
states of flavor / in event 
where (tt/ti 9 
demonstrate that 
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For independent flavors, 
if){n 9 ) 7 one may easily 



®qs = E^/^'^E " f ^• s '- ( 12 ) 



The last equality in the above equation holds solely in 
the case that Poisson statistics is applicable to the inde- 
pendent flavors (i.e., the variance a 2 is equal to the mean 
(n)) e.g., in the case where the masses of the various fla- 
vors exceeds the temperature. 

If isospin symmetry is maintained by the system then 
a 2 = a 2 when /, g belong to the same isospin multi- 
plet. The associated physical picture is that fluctuations 
of isospin, or of a product quantum number involving 
isospin are brought about by fluctuations in the popula- 
tions of flavors that carry isospin. If all the members of 
an isospin multiplet have the same mass and there exist 
no chemical potentials which favor one species over an- 
other, the fluctuations of carriers with opposing values of 
J3 compensate for each other. As a result , one obtains 
the equalities, 
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i=-J 



I 3l a 2 ~ T ^ S ^ ~ E J 3< 5 *°i " °' ( 13 ) 
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i=-I 



where, the sum is restricted to lie in a given isospin mul- 
tiplet. One uses the notation of an approximate equality 
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(~) as opposed to an exact equality (=) to high-light 
the fact that such relations hold only in the case of ex- 
act isospin invariance. However, even in the case with a 
small mass split between the u and d quark, the equal- 
ities are still approximately true and will be applied in 
this spirit. In actual calculations involving quarks and 
hadrons, in the upcoming sections, the physical masses 
of all known species will be used. Based on the above, 
the following simplifying relations between the various 
covariances may be easily derived, 



a QS = a 
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for the two 



In the definitions introduced in Refs. [21U22 
coefficients: Cbs = — ^ a %sl a s an d Cqs = 3(Tq S /v s , 
one may use Eq. (|14f) to obtain the following simplifying 
relation, 



Cqs 



3 - Cbs 



(16) 



The validity of the above equation is amply demon- 
strated by Fig. n The solid red and blue data points 
are taken form the lattice calculations of Ref. [23] , where 
both these quantities were computed independently of 
each other. The hazed cyan points represent an estima- 
tion of Cqs from the Cbs points using Eq. JEJ. In the 
lattice computation, the masses of the u and d quarks are 
set equal to each other, i.e., the lattice calculation dis- 
plays exact isospin symmetry, clearly demonstrated by 
the exact coincidence of the hazed-cyan and solid-blue 
circles. It should be reiterated that Cqs may be calcu- 
lated given a Cbs using Eq. l(Trj)l. 

Similarly, one may also define another set of related 
observables, 
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We refrain from ascribing any overall normalization fac- 
tor, as in the cases of Cbs, Cqs- none is self-evident. 
Both these quantities may be estimated on the lattice, 
in models, as well as in an experiment. In all such 
cases, independent of the phase of matter involved, one 
will recover the equality imposed by isospin symmetry 
(Eq. O), 



Cqb — 



I + Csb 



(18) 



An interesting situation is afforded for flavor SU(2), i.e., 
in a theory without strangeness. While Cbs,Csb and 
Cqs are undefined, Cqb = \ (a similar situation is that 
of the canonical ensemble where S is held fixed) . This 
last value serves as an important test of any model de- 
vised to reproduce the properties of any phase of strongly 
interacting matter and will be used in the upcoming sec- 
tions where we will compare model predictions to results 
from lattice simulations. 
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FIG. 1: (Color Online) A test of the formula Cqs 
Both ratios, Cbs and Cqs, computed on the lattice in 
Ref. [2^1 . are presented as the solid red and blue points. The 
Cbs points are then used to calculate Cqs resulting in the 
hazed cyan points. 



III. OPERATOR RELATIONS 

In the previous section, we outlined various proper- 
ties and relations between the diagonal and off-diagonal 
susceptibilities of heated strongly interacting matter. In 
this section, we extend the discussion to the operator 
structure of these susceptibilities and derive general ex- 
pectations for the value of the off-diagonal flavor suscep- 
tibilities in different quasi-particle bases. The reader not 
interested in such a study can skip to the next section, 
in which comparisons of lattice results with phenomeno- 
logical models will be made. 

In numerical simulations of lattice QCD, one evaluates 
thermal expectation values of operators weighted with 
the SU(3) gauge action S g and the fermionic determinant 
det[M], 



(O) 



J VUO{Aet[M]) n fe- s s 
J VU(dct[M]) n fe~ s e 



(19) 



In the case of the quark number susceptibility, the oper- 
ator in question is 



(20) 



N qi N q . = / d A xd A yn qz {x)n q] (y), 



where n qi (x) =: ^!i(x)^f ^i(x) : and qt,qj represent 
quarks of flavor i and j. The normal ordering removes 
the leading short distance piece, which is proportional 
to the four-volume. The operator may be decomposed 
into connected and disconnected diagrams as shown in 
Fig. [3 The locations x, y of the two operator insertions 
arc indicated by the crosses in the figures. The solid 
lines represent the valence quarks, whereas the dashed 
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FIG. 2: Connected and disconnected quark number operators. 



lines represent virtual gluons or quarks as allowed by the 
Lagrangian. If the flavors of the two quarks are the same 
(i = j), corresponding to a diagonal susceptibility, con- 
tributions emerge from both types of diagrams. Whereas 
if the two flavors are different, contributions arise solely 
from the second diagram. 

In the interaction picture, the operator of Eq. i|2U[l . for 
the off-diagonal susceptibility may be expressed in the 
simplified form, 



d 3 q t / d 3 q 3 [c 
** r,s 



"t. 



(21) 



where r, s denote the spin orientations. This expression 
may then be evaluated in a basis of weakly interacting 
partons, with the effect of the interaction introduced 
perturbatively in the basis of states. Such a starting 
point for the evaluation of the off-diagonal susceptibil- 
ity is most appropriate for the case of a weakly interact- 
ing plasma of quarks and anti-quarks, e.g., in the high 
temperature limit. In this weak coupling limit, the ex- 
pectation of the susceptibility operator may be obtained 
order by order in the strong coupling constant g. The 
leading order contribution is vanishing: 
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(N qi N q] ) - (N qt )(N q] ) 
VT 



0. 



(22) 



if % ^ j as these create and annihilate different flavors. 
The first non-zero correction to the off-diagonal suscep- 
tibility at (j, = occurs at order g 6 (see Ref. |23^. 

The opposite case is that of the strong coupling limit, 
where the system is composed of bound states of quarks 
(or anti-quarks) of flavor i with anti-quarks (or quarks) 
of flavor j. The operator of Eq. (|2L)fl may no longer be 
meaningfully evaluated perturbatively starting from free 



particle states. A convenient starting point is afforded 
by the basis of bound states: in the interest of simplicity 
we focus on the specific example of the low temperature 
phase of two-flavor QCD, i.e., a system with only u, d 
quarks (anti-quarks) and gluons. At vanishing baryon 
and charge chemical potential, the system may be ef- 
fectively described in the basis of pions 7r + ,7r~,7r°. As 
derived in the appendix, in such a system, the quark 
number operators assume the simplified forms, 



N u -- 
Nd 



(23) 



In the above, N^+,N V - is the number operator for 
7r + , 7r~ . There exists a subtle difference in the meaning of 
N q (where q stands for cither quark) and N„: the quark 
number operators are meant to indicate the net amount 
of a certain quark flavor i.e., N u is the difference between 
the number of u-quarks and u anti-quarks, whereas, the 
pion number operators N w merely indicate the number 
of pions of a certain flavor and not the difference between 
the 7r + and n~ populations. 

Within the effective form of the quark number oper- 
ators, given by Eq. 112.' '>[| . the off-diagonal quark number 
covariance in two-flavor QCD may be constructed as 



'ud 



= (N u N d ) - (N u )(N d ) 

= ((A r 7r + — iV 7r -)(iV 7r - — A r 7r +)) 



(24) 



In the case of a dilute pion gas in the grand canonical 
ensemble, the various flavors may be considered to be 
uncorrelated. As a result, 



(N V+ }(N„- 



(25) 



Incorporating the above approximation in the expression 
for a^ d , leads to the simplified form for the off-diagonal 
covariance in a dilute pion gas at low temperature and 
vanishing chemical potentials (we also denote such contri- 
butions as a^ d (M), where M denotes contributions solely 
from the mesonic sector; at a higher temperature, this 
will also include contributions from more massive mesons 
e.g., p,uj etc.), 



< d (M). 



(26) 



As the variance of cither pion species is always positive, 
we obtain the general result that at low temperature and 
vanishing chemical potentials, the off-diagonal suscepti- 
bility Xud = ^nd/V i s negative. This prediction has been 
verified in lattice calculations of \ud in Ref. [24[ ■ 

As the temperature of the system is raised, the pion 
populations (and populations of heavier mesons) as well 
as the fluctuations in the populations will increase, lead- 
ing to a drop in Xud- This trend will continue until sub- 
stantial baryon populations appear. The expression for 
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o\ d in the baryon sector is quite different from that in 
the meson sector, owing to the fact that there are no va- 
lence anti-quarks in a baryon, nor any valence quarks in 
an anti-baryon. Hence, one obtains a^ d in the baryon 
sector as 



= ((3N UUU + 2N uud + N udd 
SN mn -2N mi -N add -) 
x (3N ddd + 2N udd + N uud 

■ N udd 



~%N ddd — 2N udd 



- ((3N U 



2N,, 



-Wuuu - 2N„ 



x ((3N ddd + 2N udd + N uud 
~3N ddd ~2N add ~N aud )) 

= 2<7 lud + 2cr ldd + 2(7 lud + 2 



dd ' 



(27) 



As a result, the baryon contribution to the off-diagonal 
susceptibility is always positive. As the temperature of 
the system is raised, the contributions from mesons and 
baryons begin to compensate each other. In a weakly 
interacting hadron gas the two contributions are addi- 
tive: a\ d = a^ d (M) + a^ d (B). The increasing density of 
states in the baryon sector relative to the meson sector at 
higher energies, as well as the larger prefactors involved 
in (J^ d (B) lead to an increasing cancellation between the 
two contributions as the temperature is raised. 

In lattice computations of the temperature dependence 
of Xud one notes an initial drop followed by a rise to 
zero at T — > T c . If the picture of a weakly interacting 
hadron gas remained valid past T = T c , a^ d would con- 
tinue to rise to larger positive values. The absence of 
such behavior is an indication of the breakdown of the 
picture of a weakly interacting hadron gas near and be- 
yond T = T c . Further comparisons between the behavior 
of the off-diagonal susceptibility as well as its derivatives 
as computed on the lattice, with expectations within the 
picture of a weakly interacting hadron gas are carried out 
in the upcoming section. 



IV. LATTICE VERSUS MODELS 

As pointed out in the previous sections, the various 
susceptibilities and their ratios may be measured on the 
lattice j2f| by evaluating the average of certain opera- 
tors over a set of configurations. The appropriate choice 
of observables and their sensitivity to composite struc- 
tures was discussed in the previous section. Presently 
we focus on the results obtained from such a calculation 
on the lattice and use it to isolate the subset of models 
which describe the emergent degrees of freedom at var- 
ious temperatures in strongly interacting matter. The 
models used will be rather empirical and will require very 
little beyond arguments based on general symmetry prin- 
ciples. In all comparisons, the focus will lie expressly on 
two regions: the region below T c where one expects a 



hadronic resonance gas to be the correct set of degrees of 
freedom and above 1.5T C where one would expect to be 
firmly in the deconfincd phase. 



A. Hadron gas to quasi-particle plasma 

In the region below T c the monotonic rise of the pres- 
sure and energy density with temperature, has come 
to be understood in the picture of a weakly interact- 
ing hadronic resonance gas |l4j . as originally introduced 
by Hagedorn. This picture remains true for the case 
of the diagonal and off-diagonal susceptibilities. The 
susceptibilities are computed, assuming the condition of 
Eqs. (HU,(ITa), i.e., 



T, f ( n f) u f d f 

VT 



(28) 



where, (nt},Uf,df are the thermal average, upness and 
downness of a given hadron species /. The sum is, in 
principle, over all hadrons, but is usually truncated at 
an appropriately chosen upper mass limit. As a compar- 
ison, we plot in Fig. El, the coefficients Cbs and Cqs as 
obtained from a hadronic resonance gas spectrum, trun- 
cated at the mass of the f2 _ . As in the case of the compu- 
tation on the lattice, the plots correspond to all chemical 
potentials vanishing. One notes, that the hadron reso- 
nance gas provides a good description of the behavior of 
the ratio of susceptibilities up to the point of the phase 
transition. Here the behavior of the truncated spectrum 
fails to reproduce the sharp rise in Cbs an d the corre- 
sponding sharp drop in Cqs- It should be pointed out 
that while in the lattice simulations, exact isospin sym- 
metry has been imposed, no such condition has been re- 
quired of the hadron spectrum: the masses are taken 
directly form the particle data book. 

The results from the lattice simplify in the high tem- 
perature phase, where a general statement regarding sus- 
ceptibilities may be made: off-diagonal flavor susceptibil- 
ities are vanishing compared to the diagonal susceptibil- 
ities 0i IHl ; susceptibilities inclusive of strangeness are 
smaller compared to those which involve lighter flavors. 
This may be stated as, 



Xv 



Xds < Xud < Xs < Xd = Xv 



(29) 



where the last and the first equality applies in the case 
of isospin symmetry. For simulations at vanishing chem- 
ical potentials, the mean values of the conserved flavor 
charges are vanishing, as a result (B) = (S) = (Q) = 
Using the above equalities derived form the lattice, we 
may formulate the correlation ratio Cbs as 

_ 3 Cgg) = i( u + d + s)s) _ Xus + Xds + Xs _ j , 3Q v 



(S 2 



Models of the deconfined phase must obey the above 
constraint. The simplest model of deconfincd matter 
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where the last inequality holds in the general case. 
Csb = —1 in the case of exact SU(3)/ symmetry, i.e., 
when mass of the s-quark equals the mass of its lighter 
counterparts. Once again, it may be demonstrated that 
the simplified model of quark quasi-particles satisfies this 
requirement for the ratio of susceptibilities. Thus, from 
this standpoint, it is a viable candidate for the degrees of 
freedom of hot strongly interacting matter. Both these 
conclusions may be reduced to the single observation that 
the degrees of freedom of excited matter have to be such, 
as to display minimal strength in the covariance between 
flavors, i.e., off-diagonal flavor susceptibilities have to be 
tiny compared to the diagonal ones. 



A plasma with colored meson and diquark 
bound states? 



FIG. 3: (Color online) A comparison of the Cbs and Cqs 
calculated in a truncated hadron resonance gas at \ib — Ms = 
fiQ — OMeV compared to lattice calculations at fi — from 
Ref. (2^1 . The two hazed bands for Cbs an d Cqs for the 
hadron gas plots reflect the uncertainty in the actual value of 
the phase transition temperature T c ., which is assumed to lie 
in the range T c = 170 ± lOMeV. 

is that of non-interacting quark, anti-quark and gluon 
quasi-particles. As has been demonstrated in [2j, in 
such a situation, off-diagonal susceptibilities are identi- 
cally zero as the degrees of freedom may carry only a sin- 
gle flavor. To wit, using the notation n/ = (n/) + (nf), 
for the total number of quarks and anti-quarks of flavor 
/. 

Xus + Xds + Xs n u (l x 0) + n d {l x 0) + n s {l 2 ) 
Xs n s (l 2 ) 

= I- (31) 

The zero entries in the above equation indicate that the 
up and down flavors carry no strangeness. We also work 
in the limit described in Sect. II, where the masses of 
the quasi-particles are large enough for classical Poisson 
statistics to apply. Thus a model of quark quasiparticles 
presents a Cbs which is in agreement with that derived 
form the lattice. It should be pointed out that the nature 
of the gluon sector is irrelevant in this test. The gluons 
carry no conserved flavor and are thus oblivious to any 
such constraints. As a result, such comparisons yield no 
clues to the structure of the gluon sector. 

The next independent set of ratios of susceptibilities is 
that involving the covariances of Eq. 1151 Expressions for 
Csb may be expressed as above, 

q _ _2 Xus + Xds + Xs 

Xu + Xd + Xs + ^Xus + 2Xds + ^Xud 

= -3 — > -1, (32) 

Xu + Xd + Xs 



The picture of quasi-particlc quarks is a rather simple 
solution to the constraint of Eq. (|30|l . Indeed, such a pic- 
ture of a quasiparticlc plasma has been the subject of nu- 
merous studies in weak coupling expansions [l(| d, |26j . 
Such approximations do indeed obtain a similar behavior 
as a function of temperature as compared to the lattice 
for both the off-diagonal and the diagonal susceptibil- 
ity. Recently, an alternate picture of the QGP has been 
proposed: one where a tower of bound states of quarks 
and gluons are present besides the quasiparticles them- 
selves 0, |2(j . The larger number of particles and larger 
scattering cross-sections that result in such a mixture is 
shown to account for the pressure observed on the lattice 
as a function of the temperature. The large cross-sections 
imply very short mean free paths and lend consistency to 
the macroscopic hydrodynamic picture henceforth use to 
describe the dynamical evolution of the matter produced 
at RHIC. 

However, in a plasma containing bound states of 
quarks, in the form of colored mesons, diquarks and 
quark-gluon bound states, the correlation between flavors 
is no longer negligible compared to the diagonal suscepti- 
bility. Consider a gas consisting solely of quark-antiquark 
bound states, built from three flavors of quarks. The pos- 
sible states are ud, du, us, su, ds, sd. The flavor singlets 
qq will be ignored as they carry no conserved flavor and 
thus do not contribute to any susceptibility or covariance. 
In such a system, the ratio Cbs must vanish, because all 
states have vanishing baryon number. Expressed via the 
flavor (co-)variances, the numerator of Cbs is given as 

((u + d + s)s) = -n uS - n sU - n ds - n sS 
+ n uS + n S u + rids + n sc j 
= 0. (33) 

Thus the inclusion of mesonic states changes ratios such 
as Cbs as they cause the off-diagonal susceptibilities to 
become non- vanishing (in this case Xus = ~Xs/2)- 

Quark gluon bound states contribute similarly as 
quark-antiquark quasiparticles, whereas gluonic bound 
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states make no contribution. States such as diquarks 
have a somewhat opposite effect, diquarks belonging to 
the flavor anti-triplet (ud, us, ds) contribute us — ds — 
s 2 = +1. thus leading to (we here assume ub = 0): 



_ 2n us + 2n ds + 2n us + 2n ds _ ^ 
2tl us -\- 2n ds 

whereas the states belonging to the flavor hcxaplet pro- 
duce a coefficient, 



Cbs — 



2n us + 2n ds + 2n us + 2n ds + 8n s 



2n,, 



2n n 



> 1.0. (35) 



Ref. 0| provides masses and degeneracies for the var- 
ious bound states. The masses of all the bound states 
exceed the temperature in this model, allowing us to use 
a Maxwell-Boltzmann (MB) distribution to calculate the 
populations of such states. Such a computation was car- 
ried out in Ref. |2l| at a temperature of T = 1.5T C and 
yielded a value Cbs — 0.62 quite different form the value 
of unity found on the lattice. 



C. Introduction of baryonic bound states 

These results have motivated the inclusion of a va- 
riety of baryonic states into the model outlined above 
[27| ■ We illustrate the effect of such additions on the cor- 
relation between flavors in the following simple model. 
In the interest of simplicity, the flavor group will be re- 
stricted to SU(2)/. Lattice results for susceptibilities and 
their derivatives with respect to baryon chemical poten- 
tial using dynamical quarks exist in this case p3 |. The 
model will consist of quark and anti-quark quasi-particles 
u, d and u, d, meson like bound states uu, dd, ud, du, di- 
quark states uu, dd, ud and their antiparticles as well as 
baryons uuu, uud, udd, ddd as well as the corresponding 
anti-baryons. We assume that there is no significant co- 
variance between these quasiparticles, which arc assumed 
to be massive enough for MB statistics to apply. We will 
compute general contributions to the off-diagonal suscep- 
tibility and its various derivatives. 

In such a situation, the off-diagonal susceptibility at 
vanishing chemical potential may be decomposed as 



Xud 



VT 



[-2< d - + 2n° ud + 4n° uad + 4<J , (36) 



where, as before, 2n x includes similar contributions from 
both particles and anti-particles at vanishing chemical 
potentials. It is also assumed that populations of higher 
excited states, e.g., the hexaplet of diquarks, as well as 
states lying in the baryon decuplet are included in the 
respective populations. In the remaining, we will deal 
with densities as opposed to the absolute numbers: 



P, 



V 



(37) 



Given the off-diagonal susceptibility in a range of tem- 
peratures, one obtains a temperature dependent relation 
between the baryonic and mcsonic densities, i.e., 



2pl d {T) = 2p° ud (T) + 4p uud (T) + 4pl dd (T) 

- T Xud (T,n = 0), (38) 

where p x (T) represent the densities of various quasipar- 
ticle species at temperature T and vanishing chemical 
potential. Unlike the conventional use of the term bary- 
onic density, p° uud and p Q udd denote the density of a certain 
type of baryon and not the difference between the baryon 
and anti-baryon densities (i.e., p^ ud , p1 dd do not denote 
the net baryon density). Adjusting the baryonic densities 
compared to the mesonic densities, one may obtain the 
requisite off-diagonal susceptibility Introducing a large 
enough baryon density one may engineer a vanishing Xud 
and as an extension a vanishing Xus- With such densi- 
ties, a, Cbs = Cqs = 1 may also be achieved by a plasma 
of bound states. 

To differentiate a plasma of quasi-particle quarks and 
anti-quarks which naturally produces a \ ud — > from a 
plasma of colored bound states with a similar property 
one needs to consider the derivatives of the off-diagonal 
susceptibility Xus- At finite baryon chemical potential, 
the susceptibility will vary, as the populations of the di- 
quarks and the baryons change under the influence of 
the baryon chemical potential (quark anti-quark bound 
states carry no baryon number and hence remain unaf- 
fected). One obtains: 



T X ud(T^) c -2p° uS (T) + {2p° ud (T)+4p uud (T) 

+ 4 P l dd (T)}cosh(uf3), (39) 

where (3 is the inverse temperature. The expression is 
valid in the regime where MB statistics may be used 
instead of the full quantum statistics. Differentiating 
Eq. with respect to p,(3, one obtains the relation 



r d 2 Xud 



T 



d(pP) 2 



2p° ud (T)+Ap" uud (T) + Apl dd (T) 
d 4 Xud 



T 



(40) 



The model at this stage is applicable to any situation 
where there are baryonic and mesonic degrees of freedom 
which carry well defined quantum numbers of upness or 
downness. The baryons and mesons may be colored or 
color singlets. The considerations outlined above are ap- 
plicable to all such models, including the hadron reso- 
nance gas model used below T c . 
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FIG. 4: (Color online) A plot of the mesonic and baryonic 
contributions to Xud(T,n = 0) and and a plot of the deriva- 
tives Td 2 Xud/d((i(3) 2 \ l _ l =o (which are equal to the baryonic 
contribution ps) and Td 4 Xud/d(pf3) 4 \^ = o (which is negative 
beyond T c — 170MeV and thus is inconsistent with a bound 
state interpretation of the data.). See text for details. 



In this way, one may divide the contributions to the 
off-diagonal susceptibility and its derivatives in terms 
of mesonic and baryonic contributions. Using the mea- 
sured susceptibility and its derivatives, these contribu- 
tions may be estimated as a function of the tempera- 
ture. In the lattice computations, results are expressed 
in units of T c ; we assume T c = 0.17 GeV for definitencss. 
These are plotted as the thick solid line (mesons) and 
the solid circles (baryons) in Fig. 0] These densities for 
PB = 2p° ud (T) + 4 P ° uud (T) + 4p° udd (T) and p M = 2p° ud (T) 
satisfy Eq. (J2HI and the first equality of Eq. (gOJ). The 
condition imposed by Eq. 1)41) ll has to be satisfied by 
the derivatives of the susceptibility in such a picture of 
bound states. The fourth derivative of the susceptibility 
has been plotted as the square symbols. Despite large 
error bars one notes that, while the baryon density or 
the second derivative of the susceptibility is consistent 
with Eq. (|40|l below the phase transition temperature, 
it becomes inconsistent with such a condition above the 
phase transition temperature. Above T c , T3 4 Xud/c?(/i/3) 4 
is actually negative; hence, no composite quasiparticle 
picture is compatible with such results. It has already 
been pointed out in Ref. [24| that the signs of the various 
derivatives of the susceptibility are consistent with the 
picture of a weakly interacting quasi-particle gas. 

In the above simple model, numerous approximations 
were made. Maxwell-Boltzmann statistics was used 
throughout, variances were replaced with the mean popu- 
lations of the various flavors, and masses of various flavors 
were assumed to be independent of chemical potential. 
Such approximations were made to clearly illustrate the 
central point of this article, that the measured values of 
the off-diagonal susceptibilities and their derivatives are 



inconsistent with a picture of a composite quasiparticle 
plasma. On the other hand, the susceptibility in a weakly 
interacting model of quark quasi-particles, computed in 
the hard-thermal loop approximation has been shown to 
be consistent with the diagonal, and off- diag onal suscep- 
tibilities derived from lattice simulations [23j . It has also 
been pointed out that the signs of the derivatives of the 
susceptibilities are consistent with the quasi-particle pic- 
ture |2J|. A computation of the absolute values of the 
derivatives in such a picture is currently underway. 



V. EXPERIMENTAL OBSERVABLES 

In the previous sections, a theoretical study of var- 
ious diagonal and off-diagonal susceptibilities has been 
carried out and their relations with the degrees of free- 
dom in heated strongly interacting matter has been elu- 
cidated. In the present section, our focus will lie on the 
possible measurement of such correlations in heavy-ion 
experiments. Our considerations will be restricted to the 
measurement of the ratio Cbs which in the view of the 
authors is the most favorable from an experimental point 
of view. 

It is believed that thermalized, strongly interacting 
and deconfmed matter is transiently produced in central 
heavy ion collisions at RHIC. If one divides the whole 
system into small rapidity bins, then the fluctuations of 
conserved charges within a given rapidity bin are con- 
trolled by the degrees of freedom prevalent at the tem- 
peratures achieved. As the system expands and cools, it 
reconverts to a hadronic gas prior to freeze-out. If the 
transition to the confined phase is sudden, as in the case 
of a continuous transition and the longitudinal expansion 
is sufficiently large, then the net charge in the rapidity 
bin, set in the deconfined phase is maintained through 
the hadronic phase up to freezeout. Such fluctuations 
may then be measured event by event. The two major 
hurdles in the survival of such fluctuations through the 
hadronic phase are the contamination by hadronic fluctu- 
ations and measurability of the particles sensitive to the 
partonic fluctuations. The first issue is off lesser impor- 
tance for an observable such as Cbs a s the lightest carri- 
ers of strangeness are the kaons which are much heavier 
than the temperatures reached in RHIC collisions in the 
hadronic phase. They are produced in far fewer number 
than pions and hence do not manage to diffuse through 
multiple rapidity bins in the short time available in the 
hadronic phase. The measurement of baryon number is 
the primary problem in the estimation of Cbs- 

The detector most suitable to measurements of bulk 
fluctuations at RHIC is the STAR detector. In the 
measurement of baryon-strangeness correlations, the de- 
tector has to accurately assess the baryon number and 
strangeness in a given rapidity bin in each event. As 
the STAR detector is blind to stable uncharged particles 
it cannot measure the neutron and antineutron popula- 
tions. As a result, a measurement of ass may become 



rather difficult. Based on the discussion of Sect. II, we 
present the following recourse. A new quantum number 
is constructed, 

M = B + 2I 3 , (41) 

and fluctuations of M with respect to S are studied. In 
theoretical calculations of the the quantity ums as ou t- 
lined in Sects. II and IV, one notes that the assumption 
of isospin symmetry (see Eq. I|13f) ) reduces the covariancc 
a ms to simply ass, i-e., 

(tms = {{B + 2h)S)-{B + 2h){S) 

= a BS + 2ai 3 s ^ (Jbs- (42) 

As a result, in all theoretical models with isospin sym- 
metry Cms = C B s- 

In the experimental determination, M has the advan- 
tage that it is vanishing for all particles that do not carry 
charge or strangeness, thus M = for neutrons, antineu- 
trons, neutral pions, etc.. The experimental measure is 
thus 

En^ (n) S (n) -(£nM ( " ) )(En<7 ( " ) ) 

^ms — —o ; -9 ■ 

EJ^) 2 -(E„^) 

(43) 

In the above equation M^ n > and S^ n > are the total M and 
total strangeness within the given rapidity bin in event 
(n). One may not make the simplification of counting 
the product quantum number MS for individual flavors 
as in Eq. Ijllfl as the fluctuations in M and S are set in 
the partonic phase and the final hadrons are the result 
of decay from the deconfined phase. Hence, the differ- 
ent flavors are no longer uncorrelated as assumed in the 
derivation of Eq. . 

The presence of I3 in the observable, introduces a new 
problem in the experimental measurement. The lightest 
carriers of M are the charged pions which are numerous 
in the hadronic phase and may lead to contamination 
of the conserved charge in the chosen rapidity bin from 
neighboring bins. However, as the central rapidity bins 
at RHIC arc practically charge neutral, the possibility of 
contamination by charged pion fluctuations is greatly re- 
duced. One may divide the measured correlation between 
I3 and S into a genuine correlation and a contamination, 

VhS=°f? 8 + o%#. (44) 

As the fluctuations that result in a-j°g t are driven by pions 
in the hadronic phase, a sum over a relatively large num- 
ber of events will lead to this quantity becoming rather 
small compared to the actual correlation ctJ^j if viola- 
tions of isospin symmetry are negligible. This condition 
should hold for the produced hadronic phase over a range 
of rapidities at RHIC. This effect is illustrated in Fig. [3 
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FIG. 5: (Color online) A comparison of the two related ra- 
tios of variances Cbs, Cms as a function of the acceptance in 
rapidity from — \y max \ to |y mM |. 



where both Cbs and Cms are calculated from model 
simulations using the HIJING code psj . 

In Fig.0the correlations Cbs and Cms are estimated 
in a central Au — Au event at y/s = 200 ^4GeV. The 
acceptance in rapidity ranges from -\y max \ to \y max \, 
hence a larger y max indicates a larger acceptance. In an 
effort to further mimic the experimental acceptance, Kl 
mesons are ignored, and K$ mesons are identified either 
as a K° or a K° with 50% probability for either case. The 
results are presented as a function of y ma x- One notes 
that over the range of y m ax the two correlations Cbs and 
Cms arc rather similar. The increased fluctuations in 
isospin are the cause of the slightly larger value of Cms as 
compared to Cbs- This bodes well for the measurement 
of Cbs m RHIC experiments via a measurement of the 
quantity Cms, over a range of rapidity intervals. 

VI. CONCLUSIONS 

Both heavy ion collisions and lattice simulations of 
QCD at finite temperature present components in the 
study of heated strongly interacting matter. In this ar- 
ticle we have demonstrated that off-diagonal susceptibil- 
ities and ratios of susceptibilities have the ability to dis- 
cern the prevalent flavor carrying degrees of freedom in 
heated strongly interacting matter. The latter quantity 
may be measured on the lattice as well as in heavy-ion 
collisions, under the assumption that event-by-event fluc- 
tuations in a heavy-ion collision are set in the deconfined 
phase and are maintained through the hadronic phase. 

In this paper, the behavior of a number of observ- 
ables based on the ratio of susceptibilities Cbs, Cqs, etc. 
was explored both in the confined as well as the decon- 
fined phase. Under the assumption of isospin symmetry, 
simplifying relations between such observables were de- 
rived, which reveal the interdependence of such ratios, 
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e.g., Eqs. (|16I18|) . Results of the computations of such 
quantities in a hadronic resonance gas as well as in a non- 
interacting plasma of quarks and gluons were compared 
with calculations on the lattice (Fig. |3J). Such compar- 
isons demonstrate that the flavor carrying sector of QCD, 
is consistent with a deconfmed plasma of weakly inter- 
acting quarks and anti-quarks above T c . Below T c the 
behavior of Cbs and Cqs is consistent with that of a 
hadron resonance gas. 

The various relations between the ratios of suscepti- 
bilities relate the behavior of Cqs to that of Cbs- The 
behavior of Cbs above and below T c is caused primar- 
ily by the vanishing of the the off-diagonal susceptibil- 
ity Xus = Xds at T > T c and the negative value of its 
expectation below T c . A similar behavior is shown by 
other off-diagonal susceptibilities such as Xud- The re- 
mainder of our study focused on the behavior of the two 
flavor off-diagonal susceptibility Xud, as calculations of 
the temperature dependence of Xud as well as its vari- 
ous derivatives have been carried out in full unquenched 
lattice simulations. 

In Sects. Ill and IV, it was demonstrated that the be- 
havior of Xud as well as its various derivatives is consistent 
with that of a hadron gas below T c and a weakly inter- 
acting plasma of quarks and anti-quarks above T c . The 
behavior of the various derivatives of Xud with baryon 
chemical potential was shown to be inconsistent with a 
bound state picture above T = T c . Such an inconsis- 
tency remained even if baryonic bound state populations 
above T c were artificially enhanced to be consistent with 
the Cbs and Cqs measured on the lattice. 

Finally, we have proposed new experimental observable 
Cms related by isospin symmetry to Cbs- This observ- 
able is blind to uncharged and non-strange particles. We 
showed that it is equivalent to Cbs and may be mea- 
surable in experiments at RHIC. Estimates of Cms and 
Cbs hi HIJING simulations demonstrate the similarity 
of the two quantities over a range of rapidities at RHIC. 
This bodes well for its use as an experimental proxy for 
Cbs- Such measurements offer the possibility to directly 
probe the degrees of freedom in the deconfmed matter 
produced in high-energy heavy-ion collisions. 
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VII. APPENDIX 

In this appendix, we outline the derivation of Eq. (|23|) . 
Imagine strongly interacting matter at low temperature, 
confined in a box of volume V. The temperature is as- 
sumed to low enough for the prevalent degrees of freedom 
to be a dilute gas of pions. The state vector representing 



a pion with momentum p = yfrs ma y be expressed as 
l<> = E *f({n},{n},{m},{m},{l}) 

{n},{n},{m},{m},{l} 



x |{n},{n},{m},{m}) 




® |{1}>. (45) 

In the above equation, the vectors of integers 
{n}({n}), {m}({m}), represent the set of occupation 
numbers in different momentum states of u-quarks ( u- 
anti-quarks) and (i-quarks ( d-anti-quarks), i.e., 

{n} = {711,712,713,- • •}. (46) 

Values for rii may be or 1. The vector {1} represents the 
occupation numbers of the gluon sector and is a vector 
of integers k > 0. In Eq. g5J, |{n}, {n}, {m}, {m}) 
represents the general state vector of the quark (anti- 
quark) sector. The function ^({n}, {n}, {m}, {m}, {1}) 
represents the wave function of the one pion state with 
the constraint that the total momentum residing in this 
sector is p. Hence, 

*f({n},{n},{m},{m},{l}) (47) 
= {n},{m},{m},{l}) 

X 6\fi- E /'"•<"• + Pu^i + Pd,i m i + Pd,i™4 + P*9,diJ ' 

where, p u ^ is the momentum of the ith it-quark state with 
occupation rij. While not explicitly pointed out, the wave 
function also maintains over all color neutrality. 

Given the form of the one pion state, it is a trivial 
matter to formulate general expressions for multiple pion 
states. In this way, an effective basis of states at low 
temperature is constructed: 1 0) , 1 7r+ ) , 1 7r+ 7r+ ) , 1 7r+ 7r~ 2 ) 
etc. Interactions between these various states, is assumed 
to be small enough to be estimated in a perturbative 
formalism. The reader will note that the various states 
outlined above are orthogonal, given the orthogonality of 
the various states in the quark gluon occupation number 
basis used in Eq. (|45|l . Such n-pion states may also be 
expressed in an occupation number basis as above, e.g., 

|tt+) = |0i,0 2 ,--- ,0 Pl _i,l I>1 ,0p 1+1 ..-->. (48) 

One may now express the quark number operators as a 
matrix in the occupation number basis of pion states, 
i.e., 

N u = £ (49) 

{n},{m} 

|ni,n 2 , ■ • ■)(ni, n 2 , • ■ ■ \N u \mi, m 2 , ■ • ■){m l ,ro 2 , ■ • ■ | 
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Using the expression for the one pion state from Eq. (|45|l . 
we obtain the simple relation, 



JV u |0i,0 a ,- •• ,0 Pl _i,l Pl ,0 Pl+ i,---). = N u \tt+) 

= l<) (50) 

Similarly, the action of the upness operator on the one 
7r~ state may be computed to be 



Generalizing to the n-pion states, one obtains the gen- 
eral relations for the quark number operators in the basis 
of pions (in the limit that the pion gas is dilute i.e., the 
interactions between the different states are small ), 



N u = N„ + 
N d = N„- 



N, 



7T _ J 



(52) 



N u\n pi ) = -\tt Pi ) 



(51) 
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